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PRELIMINARY STUDY OF A WALL- STABILIZED CONSTRICTED ARC 


By Randolph A. Graves and William L. Wells 
Langley Research Center 

SUMMARY 

An iterative, implicit, finite -difference numerical technique is described which is 
suitable for obtaining solutions to the governing equations for a gas flowing in an axially 
symmetric constricted- arc heater. The method is shown to provide adequate solutions 
for three cases of simple pipe flows found in the literature, and for flow in a constricted- 
arc heater by direct comparison with experimental data. The comparison with arc-heater 
data includes static pressure, arc voltage, and wall heat flux, all as a function of axial 
location, and a radial temperature profile at one axial station. The arc-heater data were 
taken with air as the test gas at a heater inlet pressure of approximately 0.40 atm and 
two flow rates of 2.2 and 4.8 g/sec. The arc currents investigated were between 377 and 
584 amperes. 

A brief description of the experimental apparatus is included. 

INTRODUCTION 

A study of the axially symmetric constricted-arc column was undertaken in conjunc- 
tion with a spectroscopic study of the spectral dependence of radiant heat transfer from 
hot gases. In addition to the present specific application of the wall- stabilized constricted 
arc, there is currently considerable interest in the use of the arc configuration for mate- 
rials studies, plasma dynamics, and so forth. The constricted arc is a device in which a 
gas flows parallel to the axis of an electrical discharge and is confined by the cooled walls 
of the tube. The temperatures attained within the arc column are of order 10^ K and the 
processes which occur are highly complex and interrelated. Proper utilization of the con- 
stricted arc requires a knowledge of the spatial distribution of fluid mechanic, thermal, 
and electromagnetic variables. Toward this goal, many investigators have put forth sim- 
plified analytical methods which for a variety of reasons do not describe with sufficient 
accuracy the complex phenomena occurring within the hot arc column. (Refs. 1 and 2 
give an excellent review of previous work.) 

Representative of the existing numerical techniques found in the literature are those 
presented in references 2 and 3 . Both these methods have limitations which detract from 
their usefulness. Reference 2 uses an implicit method along the lines of a Crank- Nicholson 



approach which results in instabilities at arc currents greater than 300 amperes and/or 
gas flow rates greater than 3 g/sec. Reference 3 uses a fully explicit numerical tech- 
nique in which the axial step size is constrained to small values by the stability criterion. 
The development described herein is fully implicit and should be inherently more stable. 

The present work on the constricted-arc section was performed to provide a numer- 
ical tool for the description of the pertinent variables within the arc column. The numer- 
ical technique is an iterative, implicit, finite -difference method which describes the axial 
and radial fluid mechanic and thermal variables . This report will describe both the anal- 
ysis and the experimental apparatus and make comparisons between some experimental 
measurements and computed values. 


SYMBOLS 


c p 


F,f 


H 

h 

I 

k 

l 

M 


m 


N Nu 

Np r 


nodal spacing constant 

specific heat at constant pressure, J/kg-K 

axial voltage gradient, V/m 

any mathematical function 

total enthalpy, J/kg 

static enthalpy, J/kg 

electrical current, amperes 

thermal conductivity, W/m-K 

mixing length, m 

total number of nodal points 

air mass-flow rate, g/sec 


Nusselt number, 


n Nu 55 


2r wq w 
k ( T w - TJ 


Prandtl number 
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Reynolds number, 


pressure, atm (1 atm = 101.3 kN/m2) 
radiated power, W/m3 
heat flux, W/m2 


nondimensional radial distance 


radial distance, m 


temperature, K 


bulk mean temperature, T 


m = 2 C 1 UTR dR, K 


nondimensional axial velocity 


mean axial velocity, U T 


axial velocity, m/sec 




UR dR 


nondimensional radial velocity 
radial velocity, m/sec 


compressibility factor 


axial distance, m 


electrical conductivity, mhos/m 
eddy viscosity, m2/sec 
temperature, T/T m 0 
dynamic viscosity, N-sec/m2 


Planck mean absorption coefficient, 1/m 
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p 


density, kg/m3 


cr Stefan- Boltzmann constant, W/m2-K^ 

t shear stress, N/m2 

Subscripts: 
e effective 

h enthalpy 

K denotes arbitrary grid point in axial direction 

m momentum or mean 

N denotes arbitrary grid point in radial direction 

o entrance condition 

t turbulent 

w wall 

A bar over a symbol denotes a nondimensional quantity. 

APPARATUS 
Arc Heater 

The arc heater (see fig. 1) was of the wall-stabilized type, which utilizes a constant- 
diameter constrictor for stabilization and confinement of the arc and gas. The constric- 
tor was 42.4 cm long with an inside diameter of 2.37 cm. The arc attached to a hollow 
tungsten cathode at the upstream end and to multiple (32) cooled copper anodes in an 
expansion section just downstream of the constrictor exit. Resistors were in series with 
the anodes to cause the electrical current to divide evenly to each anode and thereby result 
in a very quiescent arc. The constrictor was constructed of alternate wafers of ceramic 
insulators and water-cooled heat conductors (copper). Cooling lines to these wafers were 
nonconductive so that each wafer was electrically isolated and was allowed to rise to the 
local arc potential during operation. One special wafer was equipped with a recessed 
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quartz window which provided a side view of the arc for spectrographic temperature mea- 
surements. To avoid rapid oxidation of the tungsten cathode, a very small amount (less 
than 0.2 percent of the airflow) of argon was introduced into the cathode well. 

Instrumentation 

Water and inlet air temperatures were measured with thermocouples; cooling-water 
flow rates, with turbine flow transducers; pressures, with strain- gage transducers; volt- 
age, with a voltage transductor (electrically isolated); current, with a shunt and trans- 
ductor; and the air flow rate, with a choked-flow nozzle. The data were stored on mag- 
netic tape at the rate of 400 samples per second. 

On the basis of calibrations and data- reduction techniques, the estimated accuracy 
of the measured data is as follows: 

Voltage ±3 percent 

Pressure ±2 percent 

Wall heat flux ±20 W/cm2 


ANALYSIS 
Mathematical Model 

For the purposes of the analysis, the wall- stabilized constricted arc is considered 
to be a constant-area-pipe flow apparatus; that is, flow is parallel to and symmetric about 
the tube center line (no swirl flow). No transpiration of gases at the wall is allowed; thus, 
the mass flow rate is constant. A sketch of the mathematical model of the arc heater 
showing coordinates and flow direction is given in figure 2. 

Simplifying Assumptions 

In the analysis of the wall- stabilized constricted arc, the following major assump- 
tions are made: 

(1) The flow consists of a real gas that is in local thermodynamic equilibrium. 

(2) The gas is optically transparent to all thermal radiation. 

(3) The flow is steady and axisymmetric; that is, no axial pulsations or tangential 
variations exist. 

(4) Considered as a whole, the gas is electrically neutral and the electromagnetic 
effects are negligible. 

(5) The usual velocity and thermal boundary- layer approximations are valid. (See, 
e.g., ref. 2.) 
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Governing Equations 


The geometry of the constricted arc is such that the cylindrical coordinate system 
is a natural choice for the conservation equations. Under the simplifying assumptions 
given, the Navier-Stokes equations in cylindrical coordinates for the present analysis 
reduce to: 

Continuity, 


3(pu) i 9(pvr) 
dz + r dr 


( 1 ) 


Momentum, 


ni l v 4. i JL-friL 

pn dz + 3r dz + r dr \ p dr) 


( 2 ) 


Energy, 


puP +pv |H = liJ r k 
H dz H 9r r 8 r \ c 


pc 


f + {“Pr-^“S 


0 U 


+ y&J - Q 


(3) 


where the radiated power Q is 
Q = 4p p crT 4 


In addition, Ohm's law and the integral form of the continuity equation (mass-flow equa- 
tion) will be useful in the numerical solution technique: 




yv dr 


(4) 


2tt 



pur dr = m 


(5) 


In the energy equation (eq. (3)), the Joule heating term yE z % can be simplified by using 
Ohm's law (eq. (4)) to 



( 6 ) 
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The following nondimensional variables are introduced: 


U = 


u 

u„ 


V = ^ N Re,o 


_ P-P r 
P = 


P u, 


OO 


I-f 

^0 


R = 


l TV 


Z = 


r w N Re,o 


P u 2r w 
N = IQ-O-Z 
iN Re,o fi 


xt _ ^o c p,o _ _ p 

R Pr,o jj P ~ P 

’ K o F o 


P 


Po 


k = is- 


= 


P c 


p,o 


- _ r_ 


i = i, 


N- 


Re,o 


|p o u 0 ho r w * 5 y o 


Q s Q r . wNRe ^° 

Po u o h o 


The reference conditions are calculated from 




R dR 


Upon introduction of the nondimensional quantities, the governing equations become: 
Continuity, 


£gu) l 8(p VR) .. Q 

9z + R 9R 


(7) 


Momentum, 


pjj + pY - 2l JL 
P U 8z P 3R R 9R 




( 8 ) 
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Energy, 


pU^i+pV ftt 
9z 


2 9 


N 


Pr,o 


/p k 9 h\ 1 2u 0 2 1 8 

R ^ /m , N - 1 \tT 

[ Cp &R) " N Pr>0 h Q R 9R 

K g- l w Pr,o N Pr 9R 

P 


The nondimensional form of equation (5) is 



( 9 ) 



dR = 


m 


2 ”Po u o r w 2 


( 10 ) 


Boundary Conditions 

The boundary conditions on the governing equations are 


o 

n 

> 

dH=0 

M = 0 

(at 

R = 0) 


dR 

dR 

V = 0 

o 

ii 

P 

H = H W 

(at 

R = 1) 


Turbulent Case 

When turbulent flow is considered, the transport terms are defined as follows: 


^e 


= H 


+ 


PpPtm 


k_ k , p o C P,o^m 

c p c p ^Pr.t 

e 


where 
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The eddy viscosity e m is calculated by the following relations from reference 4 : 


l = 0.4(1 - R)yL - exp 


- R V T wAy~ 


26 


e m " ^ 


dU 

dR 


(ID 


Thermophysical Properties 

The properties needed for air were assembled from several sources: 


Property 

Source 

\x, k, c p , h, Z 

Reference 5 

Mp 

Reference 6 

r 

Reference 7 


These properties as a function of pressure (or density) and temperature were input into 
the program. 


NUMERICAL SOLUTIONS 

Because of the highly nonlinear character of the governing equations, a numerical 
solution is dictated. 

In the solution scheme, the derivatives in the energy and momentum equations are 
replaced by their corresponding finite-difference form. (See the appendix.) The momen- 
tum equation contains two related unknowns, U and p; consequently, another equation 
is needed (in addition to eq. (8)) to solve for these two unknowns. Equation (10) is the 
additional relation needed, and the finite-difference form is obtained from Simpson's rule. 
(See the appendix.) In order to illustrate the solution technique, the finite-difference 
forms of equations (8) and (10) are developed as follows (see fig. 3 for nodal indexing): 
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Momentum equation (eq. ( 8 )), 


_ TT i U N,K " U N,K- 1 ^ U N+1,K ‘ (P - ^ U N,K “ U N-1,K 

p u In,k £= + p v In,k 


^6 2 + 5^ 


W + fj 


R M.j N> K + R ^-In+ 1 ,K ^ 1 n, K / 5 2 2 A 

+ / 36 2 + 6 lU 1 2 


2 8r 


U N+ 1 ,K 


5 2 2 ^ 1 n,K + R ^ 1 n- 1 ,K^ + (/3 - ^-) r m 1 n,k ^ 2 2 A + R ^In,K + R ^In+ 1 ,K 


26 , 


0 6 2 + 6 1 \5 j 5 


28 , 


U 


N,K 


6 2 2 / r M'In,K + RfX lN-l,K 


28 , 


A R P In,K / 6 2 2 A* 

1 

/ i 36 2+ 6 1 ( 6l 2 j 

U N- 1 ,K 

-J 


P K“ P K -1 


where 


6 2 “ %+l " % 


5 1 “ % " %-l 
P ~ 


»1* 


Az — z “ Z-^- ^ 


The momentum equation reduces to the desired form of 

% U N- 1 ,K + % U N,K + % U N+ 1 ,K + = % 


( 12 ) 
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where 



'N 


^ (/3 - 1 > R M-I NjK /6 2 2 \ 

| . ^N.K + ^W 


/36 2 + 6 i ^2 j 

' + 26 1 


^ V In,K 4 

R|i !n,K + R ^In+ 1 ,K ( 0 R Mn,k/6 2 2 

P 5 2 + 6 1 R L 

r n 6 2 1 + 37 

28 2 /36 2 + 6j| (g^2 


D = 

Az 


~ P K-1 + pU In,K D N,K-1 

e n' £5 

The following form results from the application of Simpson’s rule: 
Mass-flow equation (eq. (10)), 


"*1 

p UR dR 

0 


“ G 2 U 2,K + **3 U 3,K + G 4 U 4,K + R 5 U 5,K + • • • 
+ G N U N,K + %+l U N+l,K + G M-1%-1,K 


= I 


(13) 
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where 


Gn — pR 


N,K 


( 5 N + 6 N-l) 2 ( 6 N + 6 N- 1) 2 ft N + 6 N- 1 


26- 


N-l 


6 N 6 N-1 



% - pR| n K /^ 1 — S-aV^-i + In-2 «n-2 


U N-1 


2 > + 5 N + 5 N+1 


(% * 6 N - ¥ 




N 


J 


where 


6 N-2 = r N-1 ‘ %-2 


6 N-1 = % " R N-1 


6 N = %+l " % 


6 N+1 " R N+2 ' %+l 


1 = 


m 


2 »P 0 U 0 r w 2 


Equation (12) is written at each radial grid point (the grid network will be explained 
in a subsequent section), and equation (13) is added to obtain a set of simultaneous alge- 
braic equations for the unknowns Ujx and p. Additional information on this technique 
can be found in references 2,8, and 9. The matrix of coefficients of these equations 
appears as 
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The above matrix of coefficients is solved by standard techniques for the desired vari- 
ables. In a similar manner the energy equation is written in finite-difference form at 
each radial grid point, and a set of simultaneous algebraic equations which are also 
solved with standard techniques is produced. The continuity equation is solved in the 
integrated form 


1 r R d(pU) 
pR Jq dz 


dR 


rather than in finite differences to avoid the coupling of errors by the momentum and con- 
tinuity equations. (This problem was briefly discussed in ref. 2.) As a result, no stabil- 
ity problems have been experienced with the present numerical technique. 

The numerical solution scheme is fully implicit in that the solutions are iterated 
until the convergence criterion is satisfied. This procedure is shown schematically as 


13 




Nodal Grid 

The nodal spacing in the r- direction, from the wall inward to the center line, is 
taken to be a geometric progression of the form: 


r 2 = ( r l + Ar ) c r 
r 3 = ( r 2 + ir ) C r 
r i = ( r i-l + ir ) C r 

where C r is a constant usually lying between 1.0 and 1.1. The constant Ar is calcu- 
lated from 


Ar = 


(C^- 1 - l) 


c r c r 
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This particular nodal spacing allows a large number of nodal points to be spaced close to 
the tube wall, where large changes in the variables take place. Also the axial spacing can 
be uneven, as shown in figure 3 . 

Initial Profiles 

In reference 2, it was shown that the shapes of the initial profiles, at z = 0, have 
little effect on the solutions at short distances downstream; however, in the present study, 
it was found that it was necessary to give a reasonable estimate for the initial center- line 
enthalpy H^, or the effect of the initial profile would persist for some distance down the 
tube. In the present analysis, the following starting conditions are used: 

Enthalpy, 

H = H w + (H t - H w )(l.O - 3R2 + 2R3) 


Axial velocity, 

U = U t (l.O - R 2 ) 


(where is obtained from mass-flow constraints^ 


Radial velocity, 


V = 0 


RESULTS AND DISCUSSION 
Comparisons With Simple Pipe Flows 

The present numerical technique has been compared with several simplified solu- 
tions from the literature. 

Case 1: incompressible entrance flow for a circular pipe, reference 10 .- The pres- 
ent solution for laminar incompressible entrance flow and the solution of reference 10 are 
shown in figure 4 for pressure and in figure 5 for velocity. The agreement is excellent. 

Case 2: incompressible laminar-flow heat transfer in the entrance region of a cir- 
cular pipe, reference 11 .- The Nusselt numbers from the present solution and from the 
solution of reference 11 are shown in figure 6, and the solutions compare favorably. 

Case 3: compressible pipe flow with heat transfer, reference 8 .- Figure 7 shows 
the temperature profile for the present solution compared with that of reference 8. Again 
the agreement is good. An additional comparison of Nusselt numbers and ratios of wall 
temperature to mean temperature is given in table I for two nondimens ional wall tempera- 
tures. The agreement is also excellent. 
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Comparison With Arc-Heater Experimental Data 

Experimental data were taken at two different air flow rates, 2.2 and 4.8 g/sec, and 
at three different arc-current levels for each flow rate. The results of these measure- 
ments are shown in figures 8 and 9 along with the corresponding computed results. The 
parameters shown in these figures are static pressure, arc voltage, and heat flux to the 
constrictor wall, all as a function of distance from the upstream end of the arc heater. 

For an air flow rate of 2.2 g/sec, the numerical calculations closely predict the trends 
and values for all three parameters. For the air flow rate of 4.8 g/sec, the trends are 
again well predicted although the values of pressure and voltage near the downstream end 
are not in quite as close agreement with experiment as are the results for a flow rate of 
2.2 g/sec. 

In this investigation the operating current and pressure were purposely kept at low 
values in order to avoid rapid oxidation and destruction of the tungsten cathode. The val- 
ues of current were less than 600 amperes and the stagnation pressure did not exceed 
0.5 atm. According to reference 2, the electromagnetic effects are small as long as the 
current does not greatly exceed approximately 1000 amperes. At currents much higher 
than this, the electromagnetic effects may become significant and would be included in the 
governing equations. 

Likewise, the assumption of an optically thin gas is apparently valid at these pres- 
sure levels. However, at much higher pressures the validity of this assumption should be 
examined and appropriate modifications made. 

For one case (I = 377 amperes, p Q = 0.40 atm, m = 4.8 g/sec), a radial tempera- 
ture profile was measured at the axial station z = 16 cm. The temperature profile was 
determined spectroscopically by the method reported in reference 12 with the use of the 
4915 A NI line. The calculated and the measured profiles are shown in figure 10. The 
measurement and the calculations both indicate a center- line temperature of about 
12 000 K, and there is excellent agreement between the measured and the calculated 
profile over about 55 percent of the radius. The large difference between the experi- 
mental and calculated temperature values near the wall probably results from the exper- 
imental setup, primarily the recessing of the quartz window. (See ref. 12 for a discus- 
sion of the experimental procedures.) The calculated temperature profile exhibits two 
marked changes in slope corresponding to regions of rapid variation in two thermophysi- 
cal properties. The first rapid change is in the temperature range 2500 K to 3500 K, 
which corresponds to a region of rapid variation in the thermal conductivity. The second 
change is in the temperature range 5000 K to 7000 K, which corresponds to the region of 
rapid rise in electrical conductivity. Thus, what appear as oscillations in the solution are 
really thermophysical property effects. 
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CONCLUDING REMARKS 


Within the range of investigation the implicit, iterative numerical technique devel- 
oped herein for the analysis of the wall- stabilized constricted arc was shown to describe 
adequately the fluid mechanic, thermal, and electrical variables. The range investigated 
included arc currents from 377 to 584 amperes and a gas entrance pressure of approxi- 
mately 0.40 atm with gas flow rates of 2.2 and 4.8 g/sec. Predictions for operation at 
much higher pressures (several atmospheres) may require further examination of the 
assumption that the gas is optically transparent. Also, much higher currents and very 
low pressures may necessitate the inclusion of electromagnetic effects in the governing 
equations in order to maintain the degree of accuracy exhibited in this study. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 15, 1972. 
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APPENDIX 


FINITE -DIFFERENCE RELATIONS 

The following finite-difference relations are used to reduce the governing differen- 
tial equations to difference equations: 

First-order axial derivatives, 


9f_ 

3Z 


N,K 


f N,K ~ f N,K-l 
Az 


First-order radial derivatives (unequal spacing), 


_3f_ 

9R 


'N,K 


^N+1,K " ^ 


^ f N,K " f N-l,K 


/3 q 2 + 6 1 


where 



6 1 ~ % " R N-1 


6 2 - %+l " % 

Second-order radial derivatives (unequal spacing), 


JL/f 

9R\ 9R 


2 J F N + F N+1[ N+1,K ' f N,K\ 


n,k aJi ♦ Jj ’ 


+F A^- 1] 


^N+1,K " & _ 1 ^ f N,K " f N-l,K 

1 F N + F N-li 

^N,K " f N-l,K\ 

/38 2 + 

(3 2 

h 
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APPENDIX - Concluded 


Three-point unevenly spaced integration formula (Simpson's rule), 


0 

1 


o 

2 


o 

3 


f 3 f(r)clr = 


5 i + h i _ 6 i + 6 2 


36 


1 


( 5 1 


+ 6 , 



6 1 + 6 2 
36r> 


26c 


( 5 1 + 6 2) £l 5 1 + 5 2 

36o 6 0 2 
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Figure 5.- Comparison of velocity from present numerical calculations and from 
reference 10 for laminar incompressible entrance flow in a circular pipe. 
(Square symbols indicate fully developed pipe flow.) 
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Figure 7.- Comparison of temperature profile at z = 0.03 from present 
numerical calculations and that of reference 8 for compressible pipe 
flow with heat transfer. 






(a) Static pressure. 

Figure 8.- Comparison of arc-heater experimental data with numerical solutions for m = 2.2 g/sec. 
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(c) Wall heat flux. 
Figure 8.- Concluded. 



= 584 amperes 



Figure 9.- Comparison of arc-heater experimental data with numerical solutions for m = 4.8 g/sec. 










Figure 9.- Concluded. 
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